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Abstract 

The transfer matrix of the general integrable open XXZ quantum spin chain obeys 
certain functional relations at roots of unity. By exploiting these functional relations, 
we determine the Bethe Ansatz solution for the transfer matrix eigenvalues for the 
special cases that all but one of the boundary parameters are zero, and the bulk 
anisotropy parameter is wr/3, «r/5 , . . .. 



1 Introduction 



The open XXZ quantum spin chain with general integrable boundary terms £Q is a funda- 
mental integrable model with boundary, which has applications in condensed matter physics, 
statistical mechanics and string theory. Although this model remains unsolved, the special 
case of diagonal boundary terms was solved long ago [2 El Ej , and some progress on the 
more general case has been achieved recently by two different approaches. One approach, 
pursued by Cao et al. jS] is an adaptation of the generalized algebraic Bethe Ansatz [HI IZj to 
open chains. Another approach, which was developed in [IT] and which we pursue further 
here, exploits the functional relations obeyed by the transfer matrix at roots of unity. It is 
based on fusion ^2], the truncation of the fusion hierarchy at roots of unity ^3], and the 
Bazhanov-Reshetikhin ^3] solution of RSOS models. 

Both approaches lead to a Bethe Ansatz solution for the special case that the boundary 
parameters obey a certain constraint. Namely, (following the notation of the second reference 
in [TU] where 6L and a + , /3 + , 9 + denote the left and right boundary parameters, 

respectively, and N is the number of spins in the chain), 

a-+P- + a+ + P + = ±(9--0+)+7ik, (1.1) 

where k is an even integer if N is odd, and is an odd integer if N is even. This solution has 
been used to derive a nonlinear integral equation for the sine-Gordon model on an interval 
|15j . and has been generalized to other models [To] . 

Despite these successes, it would be desirable to find the solution for general values of 
the boundary parameters; i.e., when the constraint (jl.l|) is not satisfied. In the functional 
relation approach, the main difficulty lies in recasting the functional relations (which are 
known E]] for general values of the boundary parameters) as the condition that a certain 
determinant vanish. In this note we report the solution of this problem (and hence, the Bethe 
Ansatz expression for the transfer matrix eigenvalues) for the special cases that all but one 
of the boundary parameters are zero, and the bulk anisotropy has values rj = y , y , . . .. It 
may be possible to extend this analysis to more general cases. 

In Section we briefly review the construction of the transfer matrix and the functional 
relations which it satisfies at roots of unity. In Section El we present our main results; 
namely, the Bethe Ansatz solution for the transfer matrix eigenvalues when all but one of 
the boundary parameters vanish. We conclude in Section 0] with a brief discussion of these 
results. In an Appendix we briefly review the solution fOlE] for the case that the constraint 
(jl.lj) is satisfied, in order to facilitate comparison with the new cases considered here. 
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2 Transfer matrix and functional relations 



The transfer matrix t(u) of the open XXZ chain with general integrable boundary terms is 
given by [I] 



t(u) = tr K+(u) T (u) K (u) %{u) , 
where Tq(u) and Tq(u) are the monodromy matrices 

T (u) = R 0N (u) ■ ■ ■ R 01 (u) , f (u) = R 01 (u)--- R 0N (u) , 
and tr denotes trace over the "auxiliary space" 0. The R matrix is given by 

\ 

fl sinh 11 sinVi n D 

R(u) 



(2.1; 



(2.2) 



/ sinhO + r?) 

sinh it sinh 7/ 

sinh 77 sinhu 

\ sinh(M + r/) / 



(2.3) 



where 77 is the bulk anisotropy parameter; and K T (u) are 2x2 matrices whose components 
are given by |TJ H2j 



K 22( U ) 



2 (sinh a_ cosh /?_ cosh u + cosh o;_ sinh /3_ sinh u) 
2 (sinh a_ cosh /?_ cosh w — cosh a_ sinh /3_ sinh u) 



e ~ sinh 2u 



sinh 2m . 



(2.4) 



and 



(2.5) 



Ku( u ) = — 2 (sinh a + cosh /3 + cosh(u + rj) — cosh a + sinh /?+ sinh(u + 77)) 
-^22 ( M ) = —2 (sinh a + cosh /5 + cosh(u + r/) + cosh a + sinh /3 + sinh(« + i])) 
(u) = -e e + sinh 2 (w + 77) , = -e~ e+ sinh 2 (u + 77) , 

where a T , (3 T , T are the boundary parameters. 1 

In addition to the fundamental commutativity property 

[t(u),t(v)] = 0, 

the transfer matrix also has iir periodicity 

t( u + in) = t(u) , (2.7) 

1 Following |101 lllj. we use a parametrization of the boundary parameters which differs from that in 
Specifically, the matrices K^{u) are equal to those appearing in the second reference in ^H] divided 
by the factors k t , respectively. 



(2.6) 



crossing symmetry 

t(- u - rj) = t(u) , (2.8) 

and the asymptotic behavior 

e i 1 (2AT+4)+r ) (Af+2) 

t(u) ~ - cosh(#_ - B+) 1 + . . . for u — > oo . (2.9) 

For bulk anisotropy values rj = ^e_, with p = 1,2,..., the transfer matrix obeys func- 
tional relations of order p + 1 [HI CHj 

t(u)t(u + rf) . . . t(u + pr?) 

— 5(m — rj)t(u + rf)t{u + 2r/) . . . t{u + (p — l)rj) 

— 8(u)t(u + 2r])t(u + 3r]) . . .t(u + prj) 

— S(u + T))t(u)t(u + 3r])t(u + Arj) . . . t{u + prj) 

— 5(u + 2r))t(u)t(u + r))t(u + irj) . . .t(u + prj) - . . . 

— 5(u + (p - l)r/)t(u)t(u + rj) ...t(u+(p- 2)rj) 

+ ... = f(u). (2.10) 
For example, for the case p = 2, the functional relation is 

t{u)t{u + r])t(u + 2r]) - 5{u - rj)t(u + rj)- 5(u)t(u + 2r?) - 5(u + 7])t(u) = f(u) . (2.11) 
The functions S(u) and f{u) are given in terms of the boundary parameters a T , (3 T , # T by 

6{u) = SoiuMu) , f(u) = /o(«)A(u) , (2.12) 

where 

x f \ i ' v, • w i o \\ 2N sinh 2m sinh(2w + 4r^) 

<J (m) = (smhwsmh(w + 277)) — — . (2.13) 

smn(2w + rj) smh(zw + 6rj) 

5\(u) = 2 4 sinh(u + rj + a_) sinh(w + rj — a_) cosh(u + r/ + /3_) cosh(u + 77 — /?_) 

x sinh(w + 77 + aq_) sinh(w + 77 — aq_) cosh(w + 77 + cosh(w + r/ — (3 + ) , (2.14) 

and therefore, 

5(u + m)=5(u) J 5{-u - 2rf) = 8(u) . (2.15) 
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For p even, 

/ (u) = (-1) N+1 2-^ N sinh 2N ({p + l)u) , (2.16) 
/i(«) = (-l) JV+1 2 s -*( 

sinh ((p + l)a_) cosh ((p + sinh ((p + l)a + ) cosh ((p + 1)/?+) cosh 2 ((p + l)it) 

— cosh ((p + sinh ((p + cosh ((p + l)a+) sinh ((p + sinh 2 ((p + l)it) 

— (-1) N cosh ((p + 1)(0_ - 0+)) sinh 2 ((p+l)u) cosh 2 ((p+l)it)V (2.17) 

For p odd, 

/ ( w ) = (-l) 7V+1 2- 2p7V sinh 2iV ((p+l)n)t a nh 2 ((p + l)u) , (2.18) 
A(u) = -2 3 ^( 

cosh ((p + l)a_) cosh ((p + l)/9_) cosh ((p + l)ct + ) cosh ((p + 1)(3 + ) sinh 2 ((p + l)it) 

— sinh ((p + l)a_) sinh ((p + l)/3_) sinh ((p + l)ct + ) sinh ((p + cosh 2 ((p + \)u) 
+ (-1)* cosh ((p + l)(0_-0+)) sinh 2 ((p+l)u) cosh 2 ((p+l)it)V (2.19) 

Hence, /(it) satisfies 

/(« + 77) =/(«), /(-n) = /(n). (2.20) 

The commutativity property (|2.6|) implies that the eigenvectors |A) of the transfer matrix 
t(u) are independent of the spectral parameter u. Hence, the corresponding eigenvalues A (it) 
obey the same functional relations ()2.10|) . as well as the properties (12. 7|) - (|2.9|) . 

3 Bethe Ansatz solution for new special cases 

We henceforth restrict to even values of p (i.e., bulk anisotropy values rj = y , y , . . .), and 
consider the various special cases that all but one of the boundary parameters are zero. 



3.1 a- ^ 

For the case that all boundary parameters are zero except for cn_ (or, similarly, a+), we find 
that the functional relations (|2.1()|) for the transfer matrix eigenvalues can be written as 

detA4 = 0, (3.1) 
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where M. is given by the (p + 1) x (p + 1 ) matrix 

( A(u) -h(u) ... -h(-u + pq) \ 

—h(—u) A(u + prj) —h(u + pr]) ... 



M 



i —h{u + p 2 rf) ... —h(—u—p(p — l)r]) A(u + p 2 r]) J 



(whose successive rows are obtained by simultaneously shifting u i— > u + pT? and cyclically 
permuting the columns to the right) provided that there exists a function h(u) which has 
the properties 

h(u + 2m) = h (u + 2(p + l)r/) = %) , (3.3) 
% + (p + 2)77) h{-u - (p + 2)77) = 6{u), (3.4) 

P P 

Y[h{u + 2jv,)+Y[h{-u-2j V ) = f(u). (3.5) 

j=0 j=0 

To solve for h(u), we set 

h(u) = ho(u)hi(u) , (3.6) 



with 



Noting that 



W = (-l)~sinh^( M + ,) ^±M . (3.7) 

sinn(2w + 77) 



h (« + (p + 2)77) h {-u - (p + 2)77) = 6 (u), 
p p 

Y[h Q {u + 2jri) = Y[h {-u-2jri) = f Q (u) , (3.8) 

j=0 j=0 

where S (u) and fo(u) are given by ()2.13|) and ()2.16|) . respectively, we see that hi(u) must 
satisfy 

h {u + (p + 2)77) 7n(-u - (p + 2)r/) = 5 l (u) , (3.9) 
p p 

n^i(«+2j7/)+n^i(-^-2^) = (3.10) 

j=0 j=0 

Eliminating h\{—u — 2jr\) in (|3.10|) using ()3.9|) . we obtain 

p 

<«) 2 - z{u)h{u) + l[5 1 (u+ (2j - 1)77) = , (3.11) 



where 

v 

z{u) = ]Jh 1 {u + 2jr ] ). (3.12) 

3=0 

Solving the quadratic equation (|3.11|) for z(u), making use of the explicit expressions (|2.14jl 
and (|2.17|) for 8\{u) and fi(u), respectively, we obtain 



z(u) 



2~ 2(p ~ 1} cosh 2 ((p + l)u) sinh ((p + (sinh ((p + 1)«) ± sinh ((p + !)«_)) . (3.13) 



Notice that this expression for z(u) has periodicity 2r/, which is consistent with (13.12)1 and 
the assumed periodicity ()3.3|) . Corresponding solutions of (|3.12j) for h\{u) are 

2 cosh ± a_ + 77)) 

= — 4cosh u sinhw sinh(u =f a_) r-r- -r- . (3-14) 

cosh f |(w =1= a_ — 77) J 

In short, a function h(u) which satisfies ()3.3|) - ()3.5|) is given by 

h(u) = (-1) N+1 4 smh 2N (u + 77) Sinh } 2U + 27] } cosh 2 u sinh « 

, cosh (\(u + a- + 77)) 

x sinh(w-a_) — — )t7 -U- • (3.15) 

cosh [±{u — a_ — 77) J 

The structure of the matrix A4 ()3.2|) suggests that its null eigenvector has the form 
(Q(u) ,Q(u + pi]) , . . . , Q(u + p 2 r])), where Q(u) has the periodicity property 

Q(u + 2m) = Q(u) . (3.16) 

It follows that the transfer matrix eigenvalues are given by 

A{u) = h{u) — + h{-u+prj) — , 3.17 

Q(u) Q(u) 

which evidently has the form of Baxter's TQ relation. We make the Ansatz 

q( u ) = n sinn Q( m - u j)^j sinh + u i - pv^j , (3.18) 

which has the periodicity (|3.16|) as well as the crossing property 2 

Q{-u + prj) = Q(u). (3.19) 

2 Note that A(u) = A(—u + pr/) = A(—u — r/), where the first equality follows from (|3.17fl and (|3.19|) . and 
the second equality follows from the in periodicity of A(u) (which, however, is not manifest from l|3.17[) .) 
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The asymptotic behavior (|2.9|) is consistent with having M (the number of zeros Uj of Q(u)) 
given by 



M = N + p + l, (3.20) 

which we have confirmed numerically for small values of N and p. Analyticity of A(u) implies 
the Bethe Ansatz equations 

h( Uj ) _ Qfa-pn) j = l ; . . . 5 M . (3.21) 



h(-Uj + pq) Q(uj + pq) ' 

To summarize, for the special case that p is even and all boundary parameters are zero 
except for a_, the eigenvalues of the transfer matrix (|2.1|) are given by (|3.17|) . where h(u) is 
given by (|3.15j) . and Q(u) is given by (|3.18j) . (j3.2(Jj) . with zeros Uj given by (|3.21j) . 

We observe that for the special case that we are considering, the corresponding Hamilto- 
nian is not of the usual XXZ form. Indeed, t' (0) (the first derivative of the transfer matrix 
evaluated at u — 0) is proportional to a%. Hence, to obtain a nontrivial integrable Hamilto- 
nian, one must consider the second derivative of the transfer matrix. We find 



t"(0) = — 16sinh rj cosh r] sinh a 



sillll 77 

(N cosh r] + sinh ri tanh 71)0% + — : — -erf a x N ) , (3.22) 

smh a_ 



where H n n+1 is given by 



Hn ,n+l = \ + 0*0*4.1 + COsh T) . (3.23) 



3.2 



For the case that all boundary parameters are zero except for (3_ (or, similarly, /3 + ), we find 
that the functional relations (j2.10j) for the transfer matrix eigenvalues can again be written 
in the form ([3.1)1 . where now the matrix M. is given by 



M 



( A(u) -h(u) ... -h(-u-rj)\ 

—h(—u — (p+l)r]) A(u + prj) —h(u + pr]) ... 

^ -h(u + p 2 r)) ... -h(—u - (p 2 + 1)77) A(u + p 2 r)) J 



(3-24) 
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if h(u) satisfies 

h{u + 2i7i) = h(u + 2(p + 1)77) = h(u), (3.25) 

h(u + (p + 2)77) h(-u - rj) = 5{u), (3.26) 
v p 

Y[h(u + 2j v ) + l[h(-u-(2j + l)r)) = f(u). (3.27) 

j=0 j=0 

Proceeding similarly to the previous case, we now find 

h(u) = (—l) N Asinh 2N (u + r])— — \-y—^- — %- sinh 2 u cosh u ( cosh u + (— l)n sinh/JL^ . 

sinh(2w + 77) V / 

(3.28) 



The transfer matrix eigenvalues are now given by 



A{u) = h{u) — + h(-u-T]) — , (3.29) 

Q(u) Q(u) 



with 

M 



Q(u) =H sinh ( 2 



u — «,) sinh [ —(u + Uj + rj) ] , (3.30) 



which satisfies Q(u + 2«7r) = Q(w) and Q(—u — 77) = Q(u); and 

M = N + p. (3.31) 
Moreover, the Bethe Ansatz equations for the zeros Uj take the form 

h( Uj ) _ Qjuj-pn) ^ J = 1; ... ;M . ( 3.32) 
h(-Uj - 77) Q(uj +pn)' 

For this case, t'(0) = 0, and 

t"(0) = -16 cosh 77 sinh 2iV 77 K + sinh /3_ af)a%. (3.33) 

Higher derivatives yield more complicated expressions. 

3.3 e T ^ 

For the case that all boundary parameters are zero except for and 9 + (quantities of interest 
depend only on the difference 6L — 6 + ), we find that the functional relations ()2.1()j) for the 



transfer matrix eigenvalues can be written in the form (j3.1|) . where the matrix M. is given 
by 

/ A(u) -hW{-u-r)) ... -hF>(u) \ 

-hV>(u + rj) A(u + r]) -h^(-u-2ri) ... 

M = 

\ -h( 2 \-u - (p+ 1)77) ... (u + pq) A(u + pr)) J 

(whose successive rows are obtained by simultaneously shifting u 1— > u + 77 and cyclically 
permuting the columns to the right), if the functions h^ l '{u) and h^ 2 \u) satisfy 

h( k \u + m) = {u + {p + l)r]) = h {k) (u), A; = 1,2, (3.35) 

(u + v) h {2) (-u -v) = S(u) , (3.36) 

flh^iu + M + flh^i-u-jq) = f(u). (3.37) 

j=0 j=0 

We find 

hW(u) = (-1) V+- e - sinh 2 ^ + 77) Sinh ^ U + 2r >} sinh 2 2u , 

h®(u) = (-l) N e e - e + sinh 2iV (n + rj) + ^ } sinh 2 2n . (3.38) 
The transfer matrix eigenvalues are given by 

A M = fc« («)%^ + fc<»(-« - „)%+2> , (3.39) 

with, for N even, 

2M 

Q( u ) = Y[ sinh(w - tij-) , (3.40) 

which satisfies Q(w + 27r) = Q(u); and 

M=^(iV + p). (3.41) 

The Bethe Ansatz equations for the zeros Uj take the form 

hW( Uj ) _ Q(u j + r,) 

For this case, also t'(0) = 0, and 

t"(0) = -16 cosh 77 sinh 2JV 77 ^ cosh 6L cosh # + cr^ a x N + z cosh 6L sinh 9 + o\a y N 

+ i sinh 6L cosh 6* + erf cx^ — sinh 9_ sinh 6* + a\ o y N ] . (3.43) 



(3-34) 
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4 Discussion 



We have checked these solutions numerically for chains of length up to N = 6, and have 
verified that they give the complete set of 2^ eigenvalues. Hence, completeness is achieved 
more simply than in the case that the constraint (jl.lj) is satisfied 

We emphasize that, in contrast to the solution for the case that the constraint is 
satisfied, these solutions do not hold for generic values of the bulk anisotropy. Indeed, these 
solutions hold only for r/ = y , y , . . .. Also, while the Q(u) functions have periodicity m 
for the case that the constraint (jl.lj) is satisfied and for the case treated in Section EPl the 
Qiu) functions have only 2m periodicity for the cases treated in Sections 13.11 and 13.21 (See 
Eqs. (lA~T0l) . flH~UJ), (jHED and (l3~3Hj) . respectively.) 

Two key steps in our approach for solving for the function h(u) (which permits the 
recasting of the functional relations (j2.10j) as the vanishing of a determinant (j3.1jl ) are 
solving the quadratic equation (j3.11j) for z(u), and factoring the result, such as in (j3.12j) . 
For the special cases solved so far (namely, the case (jl.lj) considered in jSJ [TUJ HI] , and the 
new cases considered here), the discriminants of the corresponding quadratic equations are 
perfect squares, and the factorizations can be readily carried out. However, for general values 
of the boundary parameters, the discriminant is no longer a perfect square; and factoring 
the result becomes a formidable challenge. Perhaps elliptic functions may prove useful in 
this regard. 3 We hope to report further on this problem in the future. 
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A Appendix 

Here we briefly review the solution fHJ HTj for the case that the constraint (jl.lj) is satisfied, 
in order to facilitate comparison with the new cases considered in text. The matrix M. is 
3 An attempt along this line for the case p = 1 was considered in [Hj. 
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then given by 

I A(u) -h(-u-ri) ... -h{u) \ 

-h(u + rj) A(u + r]) -h(-u-2r)) ... 

M 



(A-i) 



y — h(— u — {p + l)rj) ... —h(u + pr]) A(u + pr]) J 
where h{u) must satisfy 

h(u + Z7r) = h (u + (p + l)r]) = h(u), (A. 2) 

h(u + rf) h(-u - rj) = 8(u) , (A. 3) 

V V 

Ylhiu + j^ + Hhi-u-jr,) = f(u). (A.4) 

3=0 j=0 

A pair of solutions is given by h(u) = h^(u) = h (u)hi (u) with h (u) given by ()3.7|) . and 
hi (u) given by 

/^(u) = (-l) JV+1 4sinh(w ± a_) cosh(w ± /3_) sinh(w ± a+) cosh(u ± (3 + ) , (A.5) 
Indeed, ho(u) satisfies 

h Q (u + 7]) h Q (-u - 77) = <5 (w), 
p p 

JJ/io^ + j^) = JJ/i (- M -jr/) = / (w), (A.6) 

j=0 j=0 

where So(u) is given by (j2.13J) . and fo(u) is given by ()2.16|) and ()2.18|) for p even and odd, 
respectively. Moreover, [1^(11) satisfies 

h^\u + T]) h ( t\-u-r ] ) = 5 1 {u), (A.7) 

where 5i(u) is given by ()2.14|) ; and 

v v 

Y[ hf\u + j V ) + Y[ hf\-u - jrj) = h(u) - (-l)^ +1 )2 1 - 2 ^sinh 2 (2(p + 1)«)) x 
i=o j=o 

x [(-1)* cosh ((p+l)(a_ + a + + (3.+(3 + )) + cosh ((p + 1)(0_ - 0+))] , (A.8) 

where fi(u) is given by (|2.17|) and (|2.19|) for p even and odd, respectively. Hence, if the 
constraint (jl.lj) is satisfied, then the RHS of (jA.8|) reduces to fi(u)] hence, all the conditions 
(|A.2j) - (jA.4|) are fulfilled. The corresponding expression for the transfer matrix eigenvalues is 
given by 
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with 

M<±) 

Q {± \u) = Yl sinh(w - uf ] ) sinh(w + uf ] + rj) , M (±) = -(N — l±k), (A.IO) 
and Bethe Ansatz equations 



h&i-uf* -77) QW^f'-ij) 
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ADDENDUM to "Bethe Ansatz derived from the functional relations of the 
open XXZ chain for new special cases" 

In [1] (to which we refer hereafter by I), we obtain Bethe Ansatz solutions for the transfer 
matrix eigenvalues of the open XXZ chain for the special cases that the bulk anisotropy 
parameter has values 

V=—r, p = 2,4,6,..., (1) 
p+1 

and one of the boundary parameters a + , is arbitrary, and the remaining bound- 

ary parameters are zero. Here we show that those results can readily be extended to the cases 
that any two of the boundary parameters {«_, a + , (3-, f3 + } are arbitrary and the remaining 
boundary parameters are either rj or in/2. (We assume that 9- = 9 + = 9.) For these cases, 
the corresponding Hamiltonians have the conventional local form (see, e.g., [2]) 

H = H n >n+ i + - sinh i] coth tanh/^Laj 1 + csch a_ sech cosh 9a\ + i sinh 9a\) 

n=l 

— coth a + tanh f3 + a z N + csch a + sech f3 + ( cosh #0"^ + i sinh #0"^) , (2) 

where i? n , n +i is given by (13.23). The corresponding energy eigenvalues are related to the 
eigenvalues A(u) of the transfer matrix t(u) (12.1) by 

E = Cl -^-A(u) +c 2 , (3) 

' " u=0 



du 



where 



1 



Cl 



16 sinh a_ cosh (3_ sinh a + cosh (3 + smh 2N 1 77 cosh r] 



smh ri + N cosh 77 . , . 

c 2 = : • (4) 

2 cosh i] 



1 a_ , a + arbitrary 

For the case that a± are arbitrary and f3± = rj, we find that 



AH 2 -4fpi («+ (2j : - 1)77) = 2~ 2p+3 cosh 2 ((p + l)u) sinh ((p+ l)u) 



x [sinh ((p + l)a_) - (-1) N sinh ((p + !)«+)] . (5) 
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The key point is that the argument of the square root is a perfect square. For definiteness, 
we henceforth restrict to even values of N . It follows that the quantity z{u) appearing in 
(13.11) is now given by (cf. (13.13)) 

z(u) = 2- 2(p - 1) cosh 2 ((j9+l)M)[sinh((p+l)M)±sinh((p+l)a_)] 

x [sinh ({p + l)u) =F sinh ((p + l)a + )\ . (6) 

Corresponding solutions of (13.12) for hi(u) are (cf. (13.14)) 

2/ cosh (\(u ± a_ + 77)) cosh (\{u =F a + + 77)) 

hi(u) = 4 cosh (u — rj) sinh(w =F «_) sinh(w ± a + ) rf- -y jf- -y 

cosh [^(u =f a_ — T])) cosh (i(-u ± a + — 7?) J 

Hence, for /i(-u) = h G (u)hi(u) we can take (cf. (13.15)) 

, / x a • i 2Ar / \ sinh(2-u + 2r/) 2 , 

mul = 4smh (w + 77) — — cosh (u — ri) 

sinh(2w + 77) 

cosh (\(u + a_ + 77)) cosh (\(u - a + + 77)) 

x smh( W -a_)smh( W + a + )- -)f) ^ J I IT' ^ 

cosh (i(-u — a_ — 77JJ cosh [±{u + a + — 77) J 

which indeed satisfies (I3.3)-(I3.5). The transfer matrix eigenvalues and Bethe Ansatz equa- 
tions are given by (13.17), (13.18), (13.21), with (cf. (13.20)) 

M = N + 2p + 1 . (9) 

2 /3- , /3+ arbitrary 

For the case that f3± are arbitrary and a± = 77, we find that 



AH 2 -4 JJ^i (u + (2j - 1)77) = ?!2~ 2p+3 sinh 2 ((p + l)u) cosh ((p + l)u) 

j=0 



x [sinh ((p + l)/3_) - sinh ((p + l)/3+)] , (10) 

and therefore 

z(u) = 2~ 2{p - 1) sinh 2 ((p+l)u) [cosh ((p+l)u)± i sinh ((p+l)/3_)] 

x [cosh((p + l)w) =F7sinh(0»+ . (11) 

Thus, we take the function /i(w) to be (cf. (13.28)) 

, / \ a ■ . 2AT/ , sinh(2-u + 27/) 2 

/i(tt) = 4smh (7/ + 77) — -sinh (u — 77) 

v 7 v #y smh(2M + 77) v " 

x (cosh 7/ + i sinh/3_) (coshw — 7 sinh/3+) , (12) 
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which indeed satisfies (I3.25)-(I3.27). The transfer matrix eigenvalues and Bethe Ansatz 
equations are given by (13.29), (13.30), (13.32), with (cf. (13.31)) 



M = N + 2p - 1 



(13) 



3 a_ , /3- arbitrary 



For the case that a_ , /3_ are arbitrary and a + = in/2, (3 + = rj, we find that 



/i(m) 2 -4 JJ^i (u+ (2j - 1)77) = 2- 2p+3 cosh 2 ((p+ l)w)sinh((p + l)u) 
3=0 



x 



sinh((p + l)a_) + (-l)ncosh ((p + (14) 



and therefore 



;(«) = 2- 2(p - 1) cosh 2 ((p+l)M)[sinh((p+l)M)±sinh((p+l)a_)] 



x 



sinh ((p + ± (— 1) n cosh ((p + 



(15) 



For /i(tt) we take 



, / \ , • 1 2AT, ,sinh(2w + 277) 
h(u) = 4 smh 27V (it + rj) ^ 



cosh(w — rj) coshw 



cosh (~(u + a_ + 77)) 

x sinh(w — a_) jf- -y (sinh w + ? cosh p_) , (16) 

cosh [±{u — a>- — rj)) 

which satisfies (I3.3)-(I3.5). The transfer matrix eigenvalues and Bethe Ansatz equations 
are given by (13.17), (13.18), (13.21), with (cf. (13.20)) 



M = N + p. 



(17) 



Similar results hold for the case that a + , f3 + are arbitrary and a_ = in/2, (5- = 77, etc. 
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